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Abstract 

We solve the Dyson-Schwinger equation for the quark propaga- 
tor in a model with singular infrared behavior for the gluon propa- 
gator. We require that the solutions, easily found in configuration 
space, be tempered distributions and thus have Fourier transforms. 
This severely limits the boundary conditions that the solutions may 
satisify. The sign of the dimensionful parameter that characterizes 
the model gluon propagator can be either positive or negative. If the 
sign is negative, we find a unique solution. It is singular at the ori- 
gin in momentum space, falls off like 1/p 2 as p 2 — ► +/ — oo, and it 
is truly nonperturbative in that it is singular in the limit that the 
gluon-quark interaction approaches zero. If the sign of the gluon 
propagator coefficient is positive, we find solutions that are, in a sense 

1 



that we exhibit, unconstrained linear combinations of advanced and 
retarded propagators. These solutions are singular at the origin in 
momentum space, fall off like 1 /p 2 asympotically, exhibit "resonant- 
like" behavior at the position of the bare mass of the quark when the 
mass is large compared to the dimensionful interaction parameter in 
the gluon propagator model, and smoothly approach a linear combi- 
nation of, free-quark, advanced and retarded two-point functions in 
the limit that the interaction approaches zero. In this sense, these so- 
lutions behave in an increasingly "particle-like" manner as the quark 
becomes heavy. The Feynman propagator and the Wightman function 
are not tempered distributions and therefore are not acceptable solu- 
tions to the Schwinger-Dyson equation in our model. On this basis 
we advance several arguments to show that the Fourier-transformable 
solutions we find are consistent with quark confinement, even though 
they have singularities on the real p 2 -axis. 

1 Introduction 

A classic approach to understanding the behavior of confined particles is to 
model and solve the Dyson-Schwinger (DS) equations for the particles' prop- 
agators. Since confinement is generally regarded as an infrared phenomenon, 
the emphasis is naturally on the infrared region of the kernels of the DS 
equations. Taking clues from studies of the infrared behavior of propagators 
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in pure Yang-Mills theory, one can adopt a vector-meson propagator model 
motivated by such studies and insert it in the kernel for the fermionic propa- 
gator equation and study issues such as fermion confinement, chiral symme- 
try breaking, the interplay between the scales for these two phenomena, and 
gauge dependence of solutions. There are two extreme views of the infrared 
behavior of the gluon propagator. One is that the singularity at q 2 = is 
much stronger than the 1/q 2 behavior of the perturbative propagator, with 
variants of 1/g 4 often proposed, and the other, in complete contrast, is that 
the propagator vanishes as q 2 — > 0. Because of the wide and rather successful 
application of the former type of behavior to bound state problems, we will 
adopt a frequently studied model of this type, proposed some time ago by 
one of the us[0, 0, for our analysis. 

Intuition for the interpretation and application of quantum field theories 
is built upon an intimate interplay between configuration space and momen- 
tum space considerations. The interaction Lagrangian and its symmetry 
properties are studied in configuration space, and space-time boundary con- 
ditions of the Green functions of the theory are crucial to their interpreta- 
tion. On the other hand, the particle spectrum and the scattering and decay 
processes contained in the theory are more intuitively assessed in momentum 
space. The particle content is revealed in the Green functions by their branch 
cuts and poles in momentum space. The complementarity of the configura- 
tion space and momentum space views is especially clear in the interchange- 
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ability of the terms short-distance and long-distance with ultraviolet and 
infrared (or hard and soft) to describe the physics of a situation. It is not 
surprising, therefore, that we assume that Green functions in (Minkowsky) 
configuration space have Fourier transforms to momentum space and vice- 
versa. Indeed, this property is at the foundation of the standard approach 
to particles and fields. 

We return to the study of fermion propagators in the infrared domain 
to see what insight can be gained by requiring that the solutions to the DS 
equation be Fourier transformable. The question is not idle, since several so- 
lutions proposed in the literature as possible models of confined behavior do 
not have Fourier transforms^, BL M. This consideration is one of the motiva- 
tions for the present work. Important collateral questions that will occupy us 
are those of the propagator behavior in the large mass limit, the asymptotic 
behavior in both timelike and spacelike directions in momentum space, and 
the behavior of the propagator in the limit that the infrared "gluon-fermion" 
interaction is turned off. We choose a simple enough model that "abelian- 
ized" Ward-Takahashi identities can be enforced at the fermion-gauge boson 
vertex and still leave us with a model whose propagator we can solve for ex- 
actly and whose Fourier transform we can evaluate. Perhaps unique to the 
present study is that we remain strictly in Minkowski space in setting up and 

solving our model equations.^ 

1 With care taken to handle the continuation to Minkowski space properly, an equivalent 
Euclidean space treatment can be given. 
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2 Defining the Model and Solving for the Prop- 
agator 

We begin by developing the model for our study of Fourier transformable 
solutions to the fermion Dyson-Schwinger equation. Our crucial ingredient 
for the DS equation is the infrared gluon propagator model 0, 0] used a 
number of times since 0, |4], [5]. As mentioned in the introduction, a number 
of studies of the gluon propagator suggest || that 

i 2 

where, in Landau Gauge, 



D{q 2 )^^ asg 2 -0 (1) 



DM) = (-9^+ q -f-)D(q 2 ). (2) 

The Fourier transform of l/(q 2 + ie) 2 does not exist 0, however, and a regu- 
larization must be prescribed to define a gluon propagator that has a Fourier 
transform. Defining D^^q 2 ) as 

where M has dimensions of (mass) 2( - A_1 \ the limit 

lim4t ) (? 2 )=^ 4 (?)(2vr)V, (4) 
where /i 2 can be positive or negative, defines the propagator for our infrared 



5 



DS equation study.Q The simple form of the confining propagator then is 

^(i)=i/iV" (5) 

in configuration space.0 

The general form of the DS equation is, in momentum space, 

l = W-m)S(p)-i J 0^^D^(k)A u (p + k,p), (6) 

where A u is the dressed vertex defined in Appendix B. Inserting our model 
propagator (4), we obtain 

l = y-m)S{p)+t/r-A ti {p,p). (7) 

We assume that A^(p + k,p) obeys the Ward identity 

Mp,p) = -^s( p ), (8) 

which is exact in an abelian gauge theory and true also in non-abelian gauge 
theory if the ghost contributions to the Ward-Takahashi identity are of order 
k and higher. The DS equation (7) then reads 

l = y- m )S(p)-iJ?r~S(p), (9) 



2 S. Blaha|| studied PP(l/(q 2 + ie) 2 ), but in a perturbative context. Pagels studied an 
alternative prescription, with a quark propagator vanishing as A — * 2, to handle the l/q 4 
singularity. This leads to a different DS Equation from ours [9]. 

3 The configuration space propagator (5) is consistent with confinement, since it clearly 
does not satisfy the cluster decomposition property [10]. 
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in our model. Taking the Fourier transform of Eq. (9), we find that the 
x-space DS equation that forms the basis for the present investigation is 

S\ x ) = (ifi - m)S(x) + z/i 2 7 • xS(x). (10) 

Phenomenologically, the dimensionful parameter \x can be expressed in terms 
of the QCD hadron dynamics scale by fitting the pion decay constant, for 
example. 

Equation (10) has the factorizable solution]^] 

Six) = e-^ 2 S (x), (11) 

where So(x) is a solution of the free DS equation. In order to avoid an expo- 
nential blow-up that invalidates the Fourier transform, we must choose an 
appropriate So(x).[] For the case /x 2 > 0, this means that the free propagator 
choice must be|]TT|| 

S (x) = -(i0 + m) (l + Ce(x )) A(x 2 ), (12) 

where C is an arbitrary constant and e(xo) A(x 2 ) obeys the free homogeneous 
DS equation. The choices C = =Fl yield advanced and retarded Green func- 
tions, respectively. A(x 2 ) obeys the inhomogeneous Klein-Gordon equation 
and has the form [11] 

A(x 2 ) = — (5{x 2 ) - —6{x 2 ) V = (13) 



4 Momentum space solutions to this model that have been offered in the literature [2,3,4] 
are not Fourier transformable. 



and 9{x 2 ) prohibits the x 2 < region where, otherwise, the full propagator 
Eq. (11) would blow up. For the /i 2 < case in Eq. (4), one needs a solution 
with 6(— x 2 ), which follows from Eq. (13) by adding the appropriate solution 
to the homogeneous equation; namelyQ 



2 , m 2 Ji(m\fx?) 



A(x 2 ) = A(x 2 ) + — ^ " ' (14) 



J i (my^x 2 ) 



my x 



with Sq(x) = + m)A(x 2 ). Note that Eq. (13) represents a tempered 
distribution, and therefore has a Fourier transform, while Eq. (14) does not. 
We emphasize that a Feynman propagator is not an acceptable choice for 
Sq(x) in Eq. (11), because the corresponding solution, S(x), does not have 
a Fourier transform. 

In summary, we have the two cases 

S(x) = e-^S (x) (15a) 

and 

S( x ) = e^Soix), (15b) 



5 This solution is the unique one in which 0(—x 2 ) appears, as required in the < 
case. That is, for this case, C — 0. 
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corresponding to the choices /i 2 > and /i 2 < 0, respectively, for the gluon 
propagator model in Eq. (4). The non-interacting Green-functions Sq(x) 
and 5o(a;)0 are manifestly regained in the fi — > limit. Avoiding the expo- 
nential blow-up as x 2 — > — oo in Eq. (15a) and as x 2 —>■ +oo in Eq. (15b) 
dictates the choices of A(x) and A(x) in Eqs. (13) and (14), as necessary 
conditions to ensure that Fourier transforms to momentum space exist. 

In our model, the Wightman function, Sw{%) =< | i/}(x)^(0) | >, 
could be identified plausibly as 

S w (x) = e-^ x \i$ + rn)W (x), (16) 



where W (x) is the free-field scalar Wightman function [11 



777 . 

W°{x) = ^K^V^m), (17) 

Any x 2 

in terms of a standard Hankel function. If Wightman functions are tempered 
distributions one can prove that free fermion asymptotic states exist flEfl ■ 
Sw{x) as defined above is not a tempered distribution, so our model is con- 
sistent with fermion confinement. The Schwinger model, which is solvable, 
illustrates such a connection between pathologies of the Wightman functions 
and confinement. The fermion Wightman functions in Coulomb gauge blow 
up exponentially in configuration space and fermion states, correspondingly, 



do not appear in the spectrum[13 



e So(x) does not have a Fourier transform. Thus (15b) is a truly nonperturbative 
solution. As we will see below, though Sq(x 2 ) is the /z — ► Q limit of (15b), the Fourier 
transform of (15b) is singular as /i — > 0. 



Let us now take up the evaluation of the Fourier transforms of Eqs. (15a) 
and (15b) and examine their behavior in momentum space. The essential cal- 
culations that must be performed are the Fourier transforms of e - ^ 2 ^ 2 A(x 2 ) 
and e + ' /i2 ' a;2 A(x 2 ); namely, choosing e~^ 2 ^ 2 A(x 2 ) for discussion, we have, 
adopting the convention fj 2 > 



and 



B{p 2 ) =-mJ d 4 xe ip - x e-^ x2 A(x 2 ){l + Ce(x )), 



(18) 



P 



'■A( P 2 ) = -ijd 



A(x 2 )(l + Ce( 



x 



= (p 2 + fi 2 p ■ d p ^B(p 2 )/m. 



(19) 



We have defined A(p 2 ) and B(p) in terms of S(p) to be 

S(p) = p'Aip 2 ) + B{p 2 ) = J d 4 xe ip - x S(x). (20) 

Details of our evaluation of the Fourier transform (18) are given in the Ap- 
pendix A, where we present a procedure that can be applied to any function 
F(x 2 ) which has a one dimensional (in x 2 ) Fourier transform. We also show 
there how to choose contours that give improved convergence for the numer- 
ical evaluation of the Fourier transform for timelike and spacelike values of 
the momentum-space argument. Writing the (7 = result for B(p 2 ) derived 
in Appendix A in the form 
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I Hi' 2 ) % 



m 



J due(u)e' ip2u+t ^tfk , (21) 



we factor out e and introduce a variable, A 2 , as follows: 



e' 1 - 2 ^ = e ^ I dr5(v — i^e 2 ^" 1 - 2 ^ 



-oo 
oo 



= e 2 ^ 2 



/ 2tt / 

— oo — oo 



<iA 2 e ^"^e 2 ^ 172 ^ 



(22) 



Substituting (22) into (21), exchanging the order of A 2 and r integration and 
evaluating the integral over v produces 

— oo 

where 



o-(A 2 ) 



e 2 ^ 



mL r dr 



g 2 M 2 l-2i(i 2 -r 

2tt 



m 



2? r /Am 



(24) 



= e^5(A 2 ) +a(A 2 )^(A 2 ) 



Thus our representation for C = is 
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B(P 2 ) _ pp \e~^ , *(A 



m 



2 

oo 2 



dA z 



, (25a) 
zr J 7r — A ZJ 
^ o ^ 

Equation (23) represents B{p 2 ) as a superposition of free propagators of mass 
A. The Fourier transform of e(x )A(x 2 , A) is obtained from that of A(x 2 , A) 
by the substitution (p 2 — A 2 ) -1 — > — 27re(]9 )<5(p 2 — A 2 ). Therefore, the Fourier 
transform of the second term in Eq. (18) is 

C[^(p )e-^^(/) + ^(/)^^/ 1 (^))e- 2 /V]. (25b) 

Equations (25a, b) show several key features of the momentum space behavior 
of the propagator that is the solution to the DS equation in our model, and 
we turn to discussion of these points in the next section. 

3 Properties of the Fermion Propagator 

3.1 The /i 2 > case 

The most obvious features of Eqs. (25a,b) are the singularities at p 2 = 0. 
As we emphasize in Sec. 4, these are not the singularitites of a Feynman 
propagator. Next we note that if /i 2 ^ 0, the A 2 integral clearly converges 
since 



h(x) 



V2 



, 1 -X 
and A 6 2M "V V 



7TX 

jr 



nm 



A3/2 



as x 



2^2 7^ 



OO 
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which is strongly convergent. Therefore the asympototic behavior of B(p 2 ) 
is 

— Til 

B(p 2 ) — ► — as p 2 ^ ±00 , (26) 

and the free-propagator ultraviolet behavior is reproduced.^ 

Next we consider the ^2 — > limit of the expressions (25a,b) for B(p 2 ). 
The singularities at the origin vanish exponentially in this limit. What hap- 
pens to the principal part integral? The asymptotic expansion for I\(x) 
shown above allows one to write the limit in the form 



a{\ 2 




2^ 



which yields 

- PP (^-^) + ™Ce( Po )5(p 2 - m 2 ), as 4 - > ( 27 ) 
m y p z — m z/ m z 

which is the Fourier transform of the free Green function (l+Ce(xo))A(x 2 ) in 
Eqs. (12) and (13). This result establishes that there is a smooth limit where 
the free momentum space Green function is the Fourier transform of the free 
configuration space Green function. This smooth limit does not obtain in our 



7 In fact, expanding Eq. (21) in powers of ^ shows that the asympototic behavior is 
given by B{p 2 )/m -> + 0\jjr) , as p 2 -> 00. 
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/i 2 < case of Eq. (4) (see Sec. (3.2) below), nor in the solutions reported 
in the literature [2, 3, 4]. B(p 2 )m is graphed for several values of /i 2 /m 2 in 
Fig. 1. The sharpening of the resonance-like behavior at p 2 = m 2 and the 
disappearance of the pole at p 2 = as ji 2 /m 2 — > is clearly shown. Thus 
in the limit as the mass of the fermion becomes large compared to the scale 
associated with the infrared behavior of the gluon propagator, the fermion 
propagator becomes more and more particle-like, in the sense that it behaves 
like l/(p 2 — m 2 ) everywhere. The pole at p 2 = is not an actual particle 
pole with the ie prescription corresponding to a time ordered product that 
insures unitarity in the perturbative expansion. This is true for any value of 
C in Eq. (12). 

A blow-up of the region near p 2 = for the 2/i 2 /m 2 = 0.2 case is shown in 
Fig. 2 to indicate just how sharp the pole is in this case where its weighting 
factor is e 1 ^ = e 5 . 

Figure 3 shows the value of (B(p 2 )/m)(p 2 — m 2 ) as a function of p 2 /m 2 
for the case 2/z 2 /m 2 = 0.2. The rapid approach to the free Green function 
behavior for large p 2 /m 2 is readily apparent. 

3.2 The Case fi 2 < - an Example of a Singular fi 2 — > 
Limit 

The free Green function A (:r 2 ), Eq. (14), and the corresponding S (x) are 
not tempered distributions and do not have Fourier transforms. Nonetheless, 
the solution (15b) to the DS equation with the vertex (8) and gluon infrared 
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propagator (4) does have a Fourier transform because the exponential factor 
e \v 2 \x 2 con trols the exp (m^\ x 2 \j divergence of Ji{m\fx 2 ) as x 2 — > — oo. 

Following the same steps as before, one arrives at Eq. (21) but with the 
opposite sign in front of [i 2 . The representation corresponding to Eq. 25a is, 

B . 9x „ r »4 1 mL m °f d\ 2 1 -2? /mA\i . . 

-p 2 =PPe^+e^- /— ^— — e ^ M— ) ■ 28 

In the limit — > ±oo, B{p 2 )/m — > ^ as in the previous case. In the limit 
fi 2 /m 2 — > 0, there is no 5(A — m) behavior, and there is no pole at p 2 = m 2 . 

m 2 

The whole expression diverges as e^ 1 in the (singular) fjr/m — > limit. 
The original free Green function A(x 2 ) is not Fourier transformable, so the 
singular nature of the /i 2 /m 2 — > limit merely reflects that fact. 

4 Discussion of Results and Conclusions 

We have reexamined a model for the infrared gluon propagator and quark- 
gluon vertex previously discussed in the literature [1-5]. We found those so- 
lutions to the quark propagator DS equation that admit Fourier transforms; 
we work directly in Minkowski space. Such solutions lend themselves to the 
study of timelike and spacelike behavior of the propagator without appeal to 
transformation to Euclidean space and continuation of the solutions found 
to the timelike region. The first solution (/i 2 > case) presented has a 
smooth limit to a combination of advanced and retarded Green functions of 
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the free Dirac equation when the interaction is turned off. In momentum 
space, the real part of the free Green function is simply PP{l/(p 2 — m 2 )), 
which is also the behavior of the full solution to the interacting model in 
the ±p 2 — > oo limit, as indicated in Fig. 3. The full propagator shows an 
interesting particle-like behavior as the mass parameter grows large com- 
pared to the infrared scale that characterizes the gluon propagator. This 
behavior is shown in Fig. 1, and it gives an explicit picture of the increas- 
ingly "free-particle-like" behavior expected as quarks become heavy. This 
is particularly true of the top quark, of course. The solution has a pole at 
p 2 = that dominates when the quark mass parameter is of the order of or 
less than the infrared scale in the gluon propagator, as shown in Fig. 1. The 
singularity at the origin is suppressed by the exp(— m 2 /2/i 2 ) factor displayed 
in Eqs. (25a,b) in the large quark mass limit. Though the model is not com- 
pletely realistic, since we do not include the ultraviolet contribution from 
the gluon propagator, it does have the interesting feature that, while the 
region near p 2 = dominates when the quark mass is small, the region near 
p 2 = m 2 dominates as the quark mass gets large. Regarding confinement, 
we note that the gluon propagator model, which is constant in configuration 
space, clearly violates the cluster decomposition property, which has been 
considered to be a sufficient condition for confinement. For the quark prop- 
agator, absence of a singularity on the real p 2 -axis is often taken to be a 
sufficient condition for confinement. In contrast, the quark propagator solu- 
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tion that we find is singular at the origin in momentum space. This singular 
behavior appears to be consistent with confinement, since the ie prescription 
necessary to define Feynman diagrams that build in the connection with uni- 
tary S-matrix elements between free outgoing and incoming colored quark 
states is not permitted by the DS Eqs. (9) or (10). Further reinforcing our 
point about confinement is the fact that the most plausible candidate for 
the two-point Wightman function in our model, Eq. (16), is not a tempered 
distribution, a condition assumed for proving the existence of free-fermion 
asympotic states. 

The second Fourier transformable solution (/i 2 < case) that we pre- 
sented has the same "free quark" asymptotic behavior in momentum space 
as the first, but the particle— like resonant-behavior at p 2 = m 2 in the large 
mass limit is missing. Furthermore, the large mass (or small interaction 
strength) limit is singular, in keeping with the fact that, while this solution 
to the full, interacting DS equation has a Fourier transform, the solution to 
the free equation does not. When taken in momentum space, the limit to 
the free Green function simply does not exist. The latter feature is shared 
by the solutions presented in [2], [3] and [4], where the (Euclidean) momen- 
tum space propagators, extended to timelike values of the argument, are not 
tempered distributions and have no Fourier transform. 

In conclusion, we have presented the tempered-distribution solutions to 
the DS equation for a simple model in 4-dimensional Minkowsky space with 
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a confining gluon propagator and a non-trivial quark gluon vertex. We offer 
the results as an interesting, instructive and useful addition to the literature 
on the dynamics of confined quarks. 
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A Appendix 

In this appendix we present details of our evaluation of the Minkowski space 
Fourier transforms. Our approach is to define the auxiliary, one dimensional 
transform, 



e~^ x A(x 2 ) = / G(u)dw, (A.l) 

— oo 

oo 

G(u) = J e-^ x2 e-^ x2 A(x 2 )dx 2 . (A.2) 



and its inverse 



It is understood that /i 2 > here. This procedure can be applied to any 
function of x 2 which has a one dimensional (in x 2 ) Fourier transform. Sub- 
stituting Eq. (13) into Eq. (A.2) and evaluating, we find 

G(uj) = J e-^^U^ 2 )--^ 2 )^^ 

— oo 

4tt 8tt J 



dx 



-4) Ji(mVx 2 ) 2 



my x 2 



1 

47T 



e -mV4(i«+^) ( A3 ) 



So the evaluation of the first term in the expression for B(p 2 ), Eq. (18), now 
involves the integral 

oo 

B(p 2 ) = ^ J due- m2/ ^ + ^ J d^xe 1 ^ 2 ^, (A.4) 

— oo 
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where we have exchanged the order of integration in writing this version of 
Eq. (18). The ^-integration is a product of Fresnel integrals that yields 

J d 4 xe tujx2+tp - x = -te(uj)^e' 1 ^ . (A.5) 

Defining v = 4/uj, the Fourier transform, B(p 2 ), can be written 

oo , 

/• 2 • 2 (1— 2ift*iQ 

o 

oo 

— / dze- z2a/i - l+a2z2) sin \tz - z/(l + a 2 z 2 )} , (A.6) 
m 2 J L J 



where a = 2p 2 /m 2 , t = p 2 /m 2 and z = z/m 2 . One obtains the Fourier 
transform of B(p 2 ) by simply changing the sign of fi 2 in (A.6). 

To evaluate B(p 2 ) for p 2 > 0, Eq. (A.6), convergence can be improved 
by choosing a contour in the first quadrant in the complex z-plane that is 
equivalent to the one shown in Eq. (A.6) by Cauchy's theorem. Writing 
z = pe td , < 9 < it/2, one can recast the representation of B(p 2 ) in the form 



n 

mB(p 2 ) = J dp[sm(xpco S 9 + 9 - ^^-)e^^ e+ ^^}, (A.7) 



p cos 9 n 

up sin i xp cos o -|- u — 

b 

where x = p 2 /m 2 , a = 2p 2 /m 2 , and D — (1 — ap sin 6 1 ) 2 + a 2 p 2 cos 2 9. The 
choice # = 7r/4 is convenient for evaluating Eq. A.7, which converges signifi- 
cantly faster than Eq. A.6 for p 2 > and thereby speeds-up the numerical 
integration. A similar formula holds for p 2 < 0, where — | < 9 < 0. 
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B Appendix 

The configuration space solution of the DS equation has the factorizable form 
S(x) = e^ x2 ^So(x) for a class of models of the vertex, independent of the gluon 
propagator's form, as we outline in this appendix. The S A (q) propagator is a 
special example where the Ward identity (8) defines the DS equation vertex 
and leads to the factorized solution (11). Let us see how this result can be 
generalized to other propagators. 



A model for the irreducible vertex that has good symmetry properties |14 
and satisfies the abelian Ward-Takahashi identity can be written in configu- 
ration space as 

T^(z; x, y) = F^y - z,x - z)S~ 1 (x, y) . (B.l) 
In Eq. (B.l) S(x,y) is the full fermion propagator, and 

F,(y -,,.-,)=/ (e-*"-» - J^ik^ + if , (B.2) 

where dF^ /dz^ = 0. One can represent the Fourier transform of in terms 
of the fermion propagator as 

T fl (p + k,p) = ^L jf 1 S-\p + ak)da + rl(p + k,p). (B.3) 

The first term on the right hand side of Eq. (B.3) is constructed to satisfy 

fcT„(p + k,p) = S-\p + k) - S-\p) , (B.4) 
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where k^Tj^(p + k,p) = 0. 

While the vertex is convenient for discussion of symmetry properties of 
a theory, the study of the fermion DS equation is generally more convenient 
when the "dressed" vertex A M is used; 

k^p + k,p) = S(p + k)T lt (p + k,p)S(p) . (B.5) 

In terms of A M (p + k, k), our model looks like 

A^(p + k,p) = jT 1 S(p + ak)da + Al(p + k,p), (B.6) 

with 

fc"A /i (p + A;,p)=5 , (p)-5 , (p + A;) (B.7) 

and /c M A^ = 0. Since A^ vanishes at k — 0, we make the approximation that 
the longitudinal term dominates in the infrared region of interest here, and 
only the first term in Eq. (B.6) is kept in what follows.^ 

The general form of the DS equation is then, in momentum space, 

l = V-m)S(p)-i J ^^D^(k)A„(p + k 7 p), (B.8) 

where D tlv {k) is the gauge boson propagator. 

With A^(p + k,p) modeled by the first term in Eq. (B.6), the Fourier 
transform of Eq. (B.8) is a pure differential equation, which we write as 

5\x) = {i0 - m)S{x) + \$f{x 2 )] S{x). (B.9) 



8 For the model of Eq. (4), does not contribute anyway, so the first term of Eq. B.6 
is the complete contribution. 
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In Eq. (B.9) we have 
d 



dx v 



f(x 2 



Xy 



di(a 2 x 2 ) + a l x l d 2 W'x 



2J i ( „2Jl\ 



(B.10) 



2xJ{x 2 ), 



with 



f(x 2 ) = -f^f(x 2 ) = - ( da\d 1 {a 2 x 2 ) + a 2 x 2 d 2 (a 2 x 2 )\. (B.ll) 
dx 2 2 Jo 1 J 

In writing Eqs. (B.8)-(10), we have expressed the covariant-gauge, gauge- 
boson propagator in the general form 



Dav{x) = d^x^g^ + d 2 (x 2 )x^x u 



(B.12) 



The general solution to the DS equation (B.9) has the remarkably simple 
factorized form 

S(x) = e if ^S (x). (B.13) 
The Green function So(x) satisfies the free equation 

(0-m)S o (x) = 8\x), (B.14) 

where it is to be understood that arbitrary solutions, Sh, to the homogeneous 
equation (ift — m)S H (x) = can always be added to solutions of Eq. (B.14). 
Our normalization of f(x 2 ) is chosen to be /(0) = 1, so the solution to Eq. 
(B.10) can be displayed as 

i 

f( x 2 ) = LJ dx' 2 J da[d 1 {a 2 x' 2 ) + a 2 x' 2 d 2 {a 2 x' 2 )], (B.15) 
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in those circumstances where the integral converges. 
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X 

Figure 1: Plot of mB(p 2 ) vs. x = p 2 /m 2 for the case C — 0, Eq. (21). 
Curves for the values a = 2/z 2 /m 2 = 0.05 (solid line), a = 0.2 (long dashed 
line) and a = 1.0 (short dashed line) are shown. The onset of the pole is 
visible for the a = 0.2 case, and it is the dominant feature in the a = 1.0 
case, where the "resonant" behavior at x = l(p 2 = m 2 ) is gone. The region 
around x = is excluded from all three plots so that all three cases fit on 
the figure (The jump across x = does not show up on this scale for the 
a = 0.05 case). 
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Figure 2: A blow-up of the p 2 = region for the a = 0.2 plot to 
display the onset of the singularity. 
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Figure 3: Plot of (x — l)B{p 2 )m for the case a = 0.2. The rapid 
approach to 1 as x = p 2 /m 2 grows, shows the l/(x — 1) asymptotic 
behavior. As seen in Fig. 1, this behavior is shared by all of 
the different a-value solutions. The same asymptotic behavior is 
obeyed by B(p 2 ) (see text). 
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